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POWER AMPLIFIER PRE-DISTORTION 



TECHNICAL FIELD 



5 The present invention relates to implementing digital pre-distortion in power 
amplifiers where memory effects occur and where parameters depend on, for 
example, the average signal power level or device temperature. 



Power amplifiers are known to add more or less distortion to the signal they 
are designed to amplify. The reason for this is that a power amplifier has a 
non-linear input-output signal characteristic. This shows up as a broadened 
spectrum around the desired amplified signal, and as an unwanted inband 

15 component of the signal. As a counter-measure to decrease the effects of 
non-linearity, it is known to pre-distort the signal at the input of the ampli- 
fier as to give an vm-distorted amplified signal at the output of the amplifier. 
This technique is called pre-distortion. Pre-distortion as implemented today 
normally uses a look-up table that is used to multiply the signal. The entry 

20 into the table is the magnitude of the signal at every time sample. 

A variation of the single look-up table approach for handling non-linearity is 
described in [1]. Here a FIR filter is used to handle both non-linearity and 
inter-symbol interference. Instead of performing the filter coefficient multipli- 
25 cation at each filter tap explicitly, look-up tables are used to perform this 
task. 

Memory effects is another problem related to power amplifiers. Memory 
efi'ects typically show up as a non- symmetrical spectrum around the carrier 
30 at the output of a power amplifier. That is, although the carrier (desired 
signal) spectrum is perfectly symmetrical, the spurious spectrum coming 
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from the distortion may be non-symmetrical with respect to the center of the 
carrier. 

The methods used to handle non-linearity do not take into account memory 
5 effects of the power amplifier. As the term "memory effects" indicates, there 
is a dependence not only on the present sample but also on previous sam- 
ples of the signal. Thus, the previously used single table approach cannot 
take care of memory effects, but can only handle non-linearity. 

10 Reference [2] suggests handling memory effects by using an envelope filter, 
which considers both the current and previous sample amplitudes in 
calculating a weighted multiplication coefficient that is intended to account 
for memory effects. 

15 Lei Ding et. al. [3] inspired by work done by Kim and Konstantinou [4] have 
derived a pre-distortion method based on what they call "Memory Polynomi- 
als" that very well model memory effects. However, this method has the 
drawback that it requires recalc\alation of the memory polynomials for each 
new input signal amplitude, which can be computationally costly, especially 

20 if many polynomials of high order are used. 

SUMMARY 

An object of the present invention is to provide a computationally efficient 
2 5 pre-distortion method based on memory polynomials. 

This object is achieved in accordance with the attached claims. 

Briefly, the present invention is based on the insight that the memory 
30 polytiomial approach can be implemented as a FIR type structxire, in which 
the "filter taps'' are replaced by look-up tables (representing sampled poly- 
nomials) triggered by the input signal amplitude. Preferably further look-up 
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tables are used to track changes in the characteristics of the power amplifier 
due to, for example, heating of semi-conductor components, 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 

The invention, together with further objects and advantages thereof, may best 
be understood by making reference to the following description taken together 
with the accompanying drawings, in which: 

Fig. 1 is a diagram illustrating the non-linear input-output signal 
1 0 characteristic of a power amplifier; 

Fig. 2 is a diagram illustrating the spectrum of the signal amplified by a 
non-linear power amplifier; 

Fig. 3 is a diagram illustrating the input-output signal characteristic of 
a power amplifier pre-distorter for removing the non-linearity in Fig. 1; 
15 Fig. 4 is a diagram illustrating the input-output signal characteristic of a 
power amplifier provided with pre-distortion; 

Fig. 5 is a diagram illustrating the spectrum of the signal amplified by a non- 
linear power amplifier with memory; 

Fig. 6 is a diagram illustrating sampling of polynomials in accordance 
2 0 with the present invention; 

Fig. 7 is a block diagram of an exemplary embodiment of a pre-distorter 
in accordance with the present invention; 

Fig. 8 is a block diagram of another exemplary embodiment of a pre- 
distorter in accordance with the present invention; 
25 Fig. 9 is a block diagram of another exemplary embodiment of a pre- 

distorter in accordance with the present invention; 

Fig. 10 is a block diagram of another exemplary embodiment of a pre- 
distorter in accordance with the present invention; 

Fig. 1 1 is a block diagram of another exemplary embodiment of a pre- 
30 distorter in accordance with the present invention; and 

Fig. 12 is a block diagram of an exemplary embodiment of a base station 
including a power amplifier provided with a pre-distorter in accordance with 
the present invention. 
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DETAILED DESCRIPTION 

In the following description the same reference designations will be used for 
the same or similar elements throughout the figures of the drawings. 

5 

Before the invention is described in detail, a brief description of the underfying 
problem will now be given with reference to Fig, 1-5. 

Fig. 1 illustrates the non-linear input-output signal characteristic of a power 
10 amplifier. At low input signal amplitudes the amplifier is almost linear, but 
at higher amplitudes it becomes more and more non-linear lontil it is satu- 
rated. This non-linearity shows up as a broadened spectrum around the 
desired amplified signal (and as an imwanted inband component of the 
signal), as illustrated in Fig. 2. As a counter-measure to decrease the effects 
15 of non-linearity, it is known to pre-distort the signal at the input of the 
amplifier to give an un-distorted amplified signal at the output of liie ampli- 
fier. This technique is called pre-distortion and is illustrated in Fig. 3. The 
input-output signal characteristic for a pre-distorted power amplifier is 
essentially linear up to saturation, as illustrated in Fig. 4. 

20 

Memory effects is another problem related to power amplifiers. Memory 
effects typically show up as a non-sjnnmetrical spectrum around the carrier 
at the output of a power amplifier, as illustrated in Fig. 5.. That is, although 
the carrier (desired signal) spectrum is perfectly sjnnmetrical, the spurious 
25 spectrum coming from the distortion may be non-symmetrical with respect 
to the center of the carrier. 

There is a theoretical way of designing a pre-distorter that takes care of all 
memory effects. This is called the Volterra series. The Volterra series is an 
30 extension to the well-known Taylor series, which can be used as a pre- 
distorter for memory- less amplifiers. The Volterra series, however, also takes 
into account time-delayed terms that may quite accurately model the pre- 
distortion, and may therefore be used to suppress the distortion spectrum. 
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However, a Volterra series quite rapidly gets large in terms of the number of 
possible terms in the expansion. For example, a polynomial of degree 5 with 
a memory depth (maximum delay) of 5 sample uriits will give rise to at least 
500 coefficients. 

5 

Since the full Volterra series can not be implemented with reasonable 
complexity, an approximation based on "Memory Polynomials" has been 
suggested in [3]. In this approximation the pre-distortion PD(n) may be 
expressed as: 

10 

PDin) = i Zak^xin - q)\xin - qf-' (1) 

where xfn-q) denotes the delayed (by q time units) samples and akq denotes 
the expansion coefficients, Unfortxmately this expression is still quite 
15 complicated, and a drawback of this prior art method is that the expression 
has to be evaluated for each new input sample x(n). However, as will be 
shown below, this expression may be rewritten into a more suitable form for 
practical implementation. The derivation essentially includes three steps: 

20 1. Separate the double sum into partial sums including only terms with 
the same delay. This gives: 
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PD(n)=X Ts^kqX(.n-q)\x(n-q)\ = 

4=1 9=0 



Separate 
delays 
into 

different 
sums 



K 



K 



+ Ti^ki^in - l)\xin - 1)|* ^ + 



K 



+ T.akQxin-Q)\x(n-Q)\ 



\k-\ 



2. Here it is noted that the delayed signals x(n-q) do not depend on the 
summation indices k. Thus the partial s\uns may be factorized into: 
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Identify the polynomials Tq(\x(n-q) \) to obtain: 



10 



PDin) = - q)Tg(\x(n - q)\) 

9=0 



(2) 
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The coefficients akq may be obtained as a least squares solution to an over- 
determined system of equations as described in [3]. 

5 It is noted in (2) that Tq(\x(n-q) \) are polynomials in the absolute value of the 
(complex) variable x(n-q). Thus, by multiplying each delayed complex sample 
x(n-q) by a polynomial in {xfrt-q) \ (which has the same delay q) and summing 
up the products for all delays q, the same end result P(n) as in [3] will be 
obtained. However, this new approach has the advantage that the polynomi- 
1.0 als Tq may be sampled at appropriate values of |x^n-q^ |, as illustrated in Fig. 
6, and stored in look-up tables. This will reduce the pre-distorter to a simple 
FIR filter structure, in which the normally constant filter coefficients are 
replaced by these look-up tables, as illustrated in Fig. 7. 

15 In the exemplary embodiment of the present invention illustrated in Fig: 7,. 
the complex input signal x(n) is forwarded to an absolute value block 10 and 
to a multiplief 12. The absolute value signal from block 10 is forwarded to a 
look-up table LUTO representing a sampled version of polynomial To. The 
corresponding (generally complex) value from look-up table LUTO is for- 

20 warded to multiplier 12, where it multiplies the input signal sample x(n). 
Input signal x(n) is also forwarded to a delay block D, where it is delayed one 
sample period for forming a delayed sample x(n-l). This delayed sample is 
processed in the same way as the non-delayed sample by an absolute value 
block 10, a multiplier 12 and a look-up table LUTl. However, look-up table 

25 LUTl now represents a sampled version of polynomial Ti instead of To. As 
illustrated in Fig. 7, further delays and look-up tables may be included (as 
indicated by the dots in the figure). Finally, the obtained products are added 
to each other in adders 14 to form the pre-distorted signal PD(n). Look-up 
tables used in accordance with the present invention make computation in 

30 real time much more efficient than the polynomial computation for each 
sample of the input signal used in [3]. The look-up tables may be updated to 
keep track of slow changes in the characteristics of the power amplifier. 
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As an example, a power amplifier with pronoimced memory effects was 
measured and then pre-distorted with a single table pre-distorter (only the 
first branch with table LUTO in Fig, 7) .The simple single table pre-distorter 
could not resolve the non-symmetiy in the measured signal spectrum at the 
5 output of the gimplifier, although it improved distortion by 5-10 dB. The 
amplifier was then pre-distorted using a double table pre-distorter (the first 
two branches with tables LUTO and LUTl in Fig. 7). This simple extension to 
a double table pre-distorter had the effect of resolving the non-sjonmetiy in 
the measured output signal. The improvement compared to a single table 
10 was in the order of 10 dB. 

A further problem associated with power amplifiers is that rapidly (but slow 
in comparison to instantaneous signal variations) changing power levels on 
the input signal will give rise to changes in power amplifier parameters. 
15 Internal heating of semi-conductor components indirectly causes this. 
Therefore, a method is sought by which we can predict these changes in a 
fast manner, and thereafter apply these adjustments. In the discussion 
above, it has been assumed that the coefficients okq in the polynomials: 



are constants. However, a more flexible approach would be to let the coeffi- 
cients depend on a slowly varying parameter ""sT representing, for exgimple, 
average input power level, transistor temperature or transistor bias (voltage, 
25 current). In such a case, equation (3) is modified to: 



20 




(3) 




(4) 



Expanding the coefficients ajtog(z) in a McLaurin series and truncating this 
30 series after M terms, one obtains: 
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m=0 



(5) 



where bhq,m are independent of z. Using this expansion, equation (4) may be 
written as: 



A:=lm=0 



(6) 



This equation may be rearranged into partial sums containing the same 
power of z\ 
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Thus, in this approximation the original polynomials Tq(\x(n-q)\) may be 
expressed as a series of pol3niomials Tqm(\x(n''q)\) multiplying different 
powers of z. The pre-distortion may in this approximation be written as; 

5 

Q 
g=0 



.7W=0 



(7) 



In this approximation a pre-distorter may be implemented by a set of look- 
up tables in each delay branch, as illustrated in Fig. 8. The expressions in 
10 the brackets Y' in equation (7) may be viewed as "filter coefficients". 



In the exemplary embodiment illustrated in Fig. 8 the delay depth is 2 (Q=2) 
and the McLaurin expansion includes 3 terms (M=3). Further terms and 
delays are possible, but these numbers have been chosen, since they are 

15 sufficiently large to illustrate the principles and sufficiently small to avoid 
cluttering of the figure. Each delay branch includes 3 look-up tables, for 
example LUTl 1, LUT12, LUT13, triggered by the same absolute signal value, 
\x(n'l)\ in this case. This embodiment also includes a block 16 for calculat- 
ing the average power z of the input signal. Block 16 may also include a 

20 smoothing filter (for example a FIR or IIR filter) to prevent abrupt changes in 
the power signal. The power value is forwarded to a multiplier 18 in each 
branch, where it is multiplied by the value from the corresponding look-up 
table LUT02, LUT12 or LUT22. The power value z is also forwarded to a 
squaring block 20 in each branch. A multiplier 22 multiplies the resulting 

25 squared power by the value from the corresponding look-up table LUT03, 
LUTl 3 or LUT23. These products are added in adders 24, and the sums are 
added to the values from look-up tables LUTOl, LUTl 1 and LUT21 in adders 
26. The resulting sums from adders 26 form the filter coefficients of the filter 
structure. From Fig. 8 and equation (7) it is noted that the embodiment of 

30 Fig. 7 and equation (2) may be considered as the special case z=0. A method 
for determining the look-up tables is described in the APPENDIX. 
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Since the parameter z is slowly varying, the average power calculated in the 
common block 16 in the embodiment of Fig. 8 may instead be calculated in 
each delay branch. Such an embodiment is illustrated in Fig. 9. Although 
5 the outputs from blocks 16 may not be identical, they are approximately the 
same due to the slowly varying average power. 

Another variation of the embodiment in Fig. 8 is illustrated in Fig. 10. This 
embodiment eliminates the squaring block 20 by rearranging the order of 
10 blocks 18, 22 and 24. It is based on the algebraic identity: 

a-^rhz + cz^ ^a + (b + cz)z 

A similar variation of the embodiment of Fig. 9 is illustrated in Fig. 11. It is 
15 appreciated that further equivalent embodiments may be obtained by other 
rearrangements of the algebraic expressions in equation (7). For example, 
instead of performing the multiplications by z and before the multiplica- 
tions by the signal and delayed signals, the multiplication order may be 
reversed. 

20 

Fig. 12 is a block diagram of an exemplaiy embodiment of a base station 
including a power amplifier provided with a pre-distorter in accordance with 
the present invention. In Fig. 12 elements that are not necessary for imder- 
standing the invention have been omitted. The baseband complex signal x(n) is 

25 forwarded to a pre-distorter 30 in accordance with the present invention. The 
pre-distorted signal is up-converted to intermediate frequency (IF) in a digital 
up-converter 32 and converted into an analog signal in a D/A converter 34, 
which in turn is up-converted to radio frequency (RF) by an analog up- 
converter 36. The RF signal is forwarded to a power amplifier 38, and the 

30 amplified signal is forwarded to an antenna. The amplified RF signal is also 
forwarded to a feedback down-conversion chain including an analog down- 
converter 40, an A/D converter 42 and a digital down-converter 44. The down- 
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converted feedback signal is forwarded to a trainer 46, which also receives the 
pre-distorted input signal for determining the look-up tables in pre-distorter 
30 in accordance with the mathematical principles described above. In this 
embodiment pre-distorter 30 may, for example, be implemented as in Fig. 8 or 
5 10, but provided with a temperature sensor 48 sensing the power amplifier 
temperature instead of a power calculation block 16. 

The pre-distorter in accordance with the present invention may be imple- 
mented as an FPGA (Field Programmable Gate Array). Another possibility is 
10 to use a micro processor or a micro/ signal processor combination and 
corresponding software. The actual computation of the look-up table entries 
may be done in an off-line manner at a slow update speed. 

If it is desirable to keep the look-up tables small, the resolution of \x(n)\ 
15 may be coarser that the resolution of x(n). 

In the description above, it has been assumed that the input signal ampli- 
tude is used as an index into the look-up tables. However, other variables 
that depend on the input signal amplitude, such as the instantaneous input 
2 0 signal power, may also be used. 

Although the present invention has been described with reference to a FIR 
filter stmcture, it is also possible to use the same principles for an IIR 
(Infinite Impulse Response) filter structure, or a combination of FIR and IIR 
25 filter structures. Thus the most general filter structure in which the inven- 
tion may be implemented is a discrete time filter structure. 



30 



It will be understood by those skilled in the art that various modifications 
and changes may be made to the present invention without departure from 
the scope thereof, which is defined by the appended claims. 
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APPENDIX 

A method to compute the tables LUTOO, LUTOl, ... is to use already stored 
tables for certain values of the parameter "2". In this APPENDIX "2" is 
assumed to represent power level. However, the method is valid also for 
other parameters, such as power amplifier (transistor) temperature. 

If we first reduce equation (7) to a memoiy-less system, this method can 
easily be explained. With the knowledge obtained from the memory-less 
system, the method can be extended to include also systems with memory. 
Thus, for a memory-less system we have: 

PDin,z) = x{n) ■ ^oi\x{n)\}+ z - T^i\xin)\j+ z"- • T2i\xin)\)\ (Al) 

The pre-distorter as expressed in equation (Al) above, is approximately valid 
over a certain range of power levels z. Now, if we compute the tables for 3 
power levels, then it should be possible to extract three new tables Tn. So, let 
us write down the three equations that will lead us to solve for Tn. Assume 
first that the table To is dedicated to 0 dBm power, then the equations will 
become: 

PD(nfi) = x(n)'[To^xin)\)] 

PDin,Zi) = xin)'[Toi\xin)\)+Zi-Tif\x(n)\)+z{^ • (l^('^)l)] (A2) 
PD(n,Z2) = x(n)\To{\x(n)\)+Z2 ■ Ti(\x(n)\)+ Z2^ •T2(\x(n}\)] 

Remember that the pre-distorters on the left-hand side of equation (A2) are 
known. By simplifying the formulation, these equations can be written in a 
slightly different way. Denoting the tables [Qo Qi Q2] at specific power levels 
[0, zi and 22] will give the following equations: 
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• Ql (N«)|) = x{n) . [To (\x(n)\)+ z, • T, z,^ . (j;c(«)|)] (A3) 
^(«)-e2(K«)|) = ^(«)-h(l^(«)|)+^2 •^i(|^(«)|)+^2^ -^ad^C")!)] 

The solution to equation {A3) is now readily at our hand as: 



'Toi\x{n)\j 




"10 0" 


-1 


"eo(k(«)|)" 






1 Zi 




QAxin)\) 






1 ^2 ^2 _ 




Q2(\x(ni) 



So concluding this part, we see that the interpolating polynomial in equation 
(Al) can directly be used to compute corrections to the basic tables com- 
puted at for example 2?=0 dBm power. The tables Tn are derived from equa- 
10 tion (A4), which in turn makes use of stored tables Qn. A new table for any 
other signal power can be computed by simply setting z^power level, and 
summing up the tables with these weights. 



Solving the inverse matrix in equation (A4) will give the following result: 

15 



'T^i\x(n)\f 




" 1 


0 0 ' 






TAx(ni) 






^22 ^23 






x(ni) 


T2\xini(} 




_^31 


C32 C33_ 






x{n)\)_ 



So, the pre-distorter in equation (Al) can also, using the inverse matrix 
elements "c" of equation (A5), explicitly be written as: 



20 



PDin,z) = x(n) 



"fio(i^(»)|)+ 

+ +C22'Z + C32 ■• z^ )• (\xin)\)- 



'■)'Q2<\xi^)\) _ 



(A6) 
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As a furtJaer refinement to the solution, exactly the same strategy can be 
used for the polynomials (or tables) for the memory system. The two systems 
can be treated separately without any coupling among them. 

5 As an example, a power amplifier is measured for input and output signals 
at 3 different power levels: 0 dB, -3 dB and —9 dB. For each power level, 
tables Qo-2 that work as pre-distorters at their respective power levels are 
computed. The objective is to compute the parameters "c", and then to use 
equation (A6) to compute a new table for power level z=-6 dB, for example. 

10 

To compute the "c"-parameters, set zi=-3 dB, z2=-9 dB in the calibration 
step, then the "c" -parameters will become: 

"1 0 0 

^21 ^22 ^23 
^31 ^32 ^33 

15 

The result of using this interpolated pre-distorter compared to use a directly 
computed pre-distorter works very well, with only a minor difference in 
spectrum level close to the carrier. 

20 It is also possible to calculate To, Ti and T2 in equation (A3) by expanding 
them into polynomials in \x(n)\ and solving for the coefficients using the 
Least Mean Squares (LMS) method. This gives the further advantage of an 
"automatic" extrapolation of the initial tables. Thus, if tables are available for 
low power levels, corresponding tables for higher power levels can be ob- 

25 tained by extrapolation according to this method. 



1 0 0 
1 (-3) (-3f 
1 (-9) (--9)^ 



1—1 



10 0 ■ 

0.444 -0.5 0.056 
0.037 -0.056 0.019 



(A7) 
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